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We study the improvement of staggered fermions using hypercubically smeared (HYP) links. We 
calculate the strange quark mass and the kaon B-parameter, Bk, in quenched QCD on a 16^ X 64 
lattice at /3 = 6.0. We find ms(MS,2 GeV) = 101.2 ± 1.3 ± 4 MeV and Ba'(MS,2 GeV) = 
0.578 ± 0.018 ± 0.042, where the first error is from statistics and fitting, and the second from using 
one-loop matching factors. The scale (1/a = 1.95GeV) is set by Mp, and rris is determined using 
the kaon mass. Comparing to quenched results obtained using unimproved staggered fermions and 
other discretizations, we argue that the size of discretization errors in Bk is substantially reduced 
by improvement. 

PACS numbers: ll.15.Ha, 12.38.Gc, 12.38.Aw 



I. INTRODUCTION 

One of the major sources of uncertainty in using preci- 
sion experimental data to constrain the standard model 
is the lack of knowledge of the matrix elements of the ef- 
fective weak Hamiltonian between hadronic states. The 
kaon bag parameter Bk , which parameterizes the matrix 
element of the AS* = 2 operator responsible for kaon- 
antikaon mixing, is one such key input for the determi- 
nation of the CKM mixing matrix. It is defined as the 
dimensionless ratio 

^ _ (j^°|57^(l-75)rf^7M(l-75)rf|g°) 

f(XO|s7^75d|0)(0|s7^75rf|^") ^ ^ 

Different approaches, including chiral perturbation the- 
ory, the large Nc expansion, QCD sum rules and lattice 
QCD, have been used to estimate Bk- The advantage 
of the lattice approach is that it is a first principle, non- 
perturbative determination. On the other hand it in- 
troduces statistical and systematic errors like those due 
to discretization and the matching of lattice and con- 
tinuum operators. To gain control over these uncertain- 
tics, different fermion discretizations — Wilson, staggered, 
domain-wall (DW) and overlap — have been used in sim- 
ulations.^ 

In this note we explore the extent to which improved 
staggered fermions can be used to reduce two of the most 
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important systematic errors, i.e., those due to discretiza- 
tion and the matching of lattice and continuum opera- 
tors. This test is carried out in the quenched approxi- 
mation to get an estimate of the size of these errors by 
comparing with existing data. Our ultimate aim, how- 
ever, is to find a method which can be used effectively 
on dynamical lattices likely to be produced in the near 
future. 

Staggered fermions are an attractive choice for the cal- 
culation of weak matrix elements because they are com- 
putationally efhcient — indeed, simulations with three dy- 
namical flavors are already possible with relatively light 
quark masses — and yet retain sufficient chiral sym- 
metry to protect operators of physical interest from mix- 
ing with others of wrong chirality. Their disadvantage is 
that they retain four "tastes" of doublers for each lattice 
field. In the continuum limit, these four tastes become 
degenerate, and one can remove the additional degrees 
of freedom by hand. For the valence quarks this proce- 
dure is explained for the calculation of Bk in Ref. 0], 
while for the sea quarks one must take the fourth-root 
of the quark determinant. At non-zero lattice spacing, 
however, quark-gluon interactions violate the taste sym- 
metry. This has three important consequences for calcu- 
lations of Bk- 

The first concerns taste symmetry violation and the 
need to take the fourth-root of the quark determinant. 
For non-zero lattice spacing, there is no proof that the 
underlying lattice action is local and lies in the same uni- 
versality class as QCD. Even though we do not face this 
problem in quenched simulations, it is relevant when ex- 
tending our calculations to dynamical simulations. Our 
justification for proceeding is empirical — accurate un- 
quenched simulations using the fourth-root of the deter- 
minant find agreement between lattice and experimental 
results 

Second, large 0{a?') discretization errors have been 



2 



TABLE I: Quark masses used in simulation and their relation 
to the strange quark mass. 



Name 


arUq 


mq/nis 


mi 


0.01 


0.192 


m2 


0.02 


0.385 


m-i 


0.03 


0.577 


mi 


0.04 


0.769 



observed in the calculation of masses and matrix ele- 
ments. Overcoming these requires the use of very small 
lattice spacings to make reliable continuum extrapola- 
tions. Lastly, many one-loop perturbative estimates of 
matching factors differ significantly from their tree-level 
value of unity, raising doubts about their accuracy Ii| . 

The purpose of this paper is to show, using and Bx 
as probes, that the latter two problems can be greatly 
alleviated by improving staggered fermions using "fat" 
links . Based on the analysis of Ref . Q , we choose a 
particular type of fattening, hypercubic (HYP) smeared 
links 7] , although we expect that other choices will work 
comparably well. Earlier calculations show that taste- 
symmetryviolations in the spectrum are substantially re- 
duced 0,13 , ^'^'i one- loop corrections to matching factors 
for four-fermion operators which were as large as 100% 
are now reduced to ~ 10% |||. The largest improve- 
ment is in renormalization constants of left-right (pen- 
guin) four-fermion operators and this can be traced 
back to the improvement in Zm — ^/Zs — l/Zp « 1 
with HYP smearing. This has a major impact on the 
extraction of as we show in section IIIII In the case 
of Bk the major impact of improvement is to reduce 
discretization errors. This is because the one-loop cor- 
rections to matching factors in this case turn out to be 
small (~ 10%) before (as well as after) improvement. 

To test the efficacy of improvement for and Bk , we 
compare our results to the JLQCD analyses with unim- 
proved staggered fermions that include detailed studies 
of both discretization and perturbative errors 

This paper is organized as follows. In Sec. [HI we 
analyze the p meson spectrum calculated on the HYP 
smeared lattices and obtain the lattice scale 1/a. In 
Sec. IIIII we present the extraction of strange quark mass 
from the pion spectrum and compare the result with 
that obtained using unimproved staggered fermions. In 
Sec. IIVI we present results for Bk calculated using the 
HYP improved staggered fermions and compare them 
with those of unimproved staggered fermions and with 
some recent data obtained using domain wall and overlap 
fermion formulations. We close with some conclusions in 
Sec. El 



II. p MESON SPECTRUM 



Wilson plaquette action at /9 = 6.0. The lattices were 
first HYP smeared using the tree-level improved param- 
eters of Ref. 0. On these HYP smeared lattices quark 
propagators are calculated using 2Z wall sources on time- 
slices 0, 16, 32 and 48 for the bare quark masses listed in 
Tabled ^ Meson correlators are calculated at all time 
slices for each set of the 2Z wall sources. Throughout 
this work we only consider mesons composed of degener- 
ate quarks. 

The lattice scale is set using the p meson mass. We 
calculated correlators for two types of p mesons: p(l), 
with spin-taste (7i(S'Ci); p(2), with spin-taste (7174 (gi 
^i^4). The correlators are fit to the standard form |13| : 

C{t) = Zi{exp[-mit]+exp[-mi{L-t)]} 

+ Z2{-iy {exp [-m2t] + exp [-m2(L - t)]}{2) 

where mi is the mass of the p meson, while m2 is the mass 
of its opposite parity partner, whose contribution has an 
alternating sign in the time direction. The partner for 
p(l) is the bi meson, with spin-taste {jjjk <8)'CjCfe), while 
that for the p{2) is the ai with spin-taste (7^75 ® £,1^5)- 

We obtain good fits to both p correlators except for the 
p{l) at the lightest quark mass. The resulting masses, 
as well as those of the parity partner ai, are given in 
Table ITU Since we want to be able to use all four quark 
masses to carry out the chiral extrapolation, we opt to 
consider only the p{2) results to determine the lattice 
scale. We do not, however, expect that the resulting scale 
would change significantly were we to use p(l) masses, 
because the p{l) and p(2) masses agree within errors for 
the three heavier quark masses. To illustrate the quality 
of the fits we show the effective mass plots for p(2) as 
a function of time in Figs. I1I4I The effective mass at 
time t = T is defined to be the value of mi obtained 
by solving Eq. Q using the correlation function on time 
slices T to r + 3. All errors are determined using single 
elimination jackknife, with the underlying fits (both in t 
and m) using uncorrelated errors, because the correlation 
matrix is determined with insufficient accuracy. 

The result of a quadratic fit to Mp(^2) versus quark 
mass, shown in Fig. [S] is 

aMp = 0.399(10)4-2. 60(59)(aTO,)-9.3(8.7)(aTO,)2. (3) 

From this we estimate the lattice scale 1/a quoted 
in Table ^1 by setting the chirally extrapolated value 
0.399(10) = aMP'^y . The change in the result- 
ing scale from extrapolating to the physical light quark 
masses {a{mu -\- rnd)/2 ^ 0.0015) rather than the chiral 

limit is smaller than our statistical errors and wc do not 

1/2 3/2 

include it. We also do not include the mq and niq 
terms from pion loops in our chiral fit, as they are 
expected to have small coefficients, and we have too few 
mass points to reliably include them. 



The statistical sample consists of an ensemble of 218 
gauge configurations of size 16^^ x 64 generated using the 



^ The details on the 2Z wall source are given in Ref. m. 
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TABLE II: Masses of p and a\ mesons, and resulting scales. 







ciMp(7i74 ® CiC4) 


aMai(7i75 ® Ci6) 


mi 




0.4244(63) 


0.5897(345) 


m-i 


0.4444(32) 


0.4466(40) 


0.6350(218) 


7713 


0.4676(32) 


0.4692(43) 


0.6387(347) 


m4 


0.4865(25) 


0.4879(32) 


0.6657(250) 


1/a 




1945(50) MeV 


2112(131) MeV 



A potential problem with our estimate of the scale 
is our use of a relatively small volume [L « 1.6 fm). 
Although we expect this is large enough to study kaon 
properties (since ttikL « 5 is larger than the range 3 — 4 
where significant effects usually set in), we are relying 
in our scale determination on results from all four quark 
masses. At the lightest quark mass M^iy — 2.7, and vol- 
ume errors may be significant. Evidence that this is the 
case comes from Ref. who have results for Mp for 
three different volumes a,t P — 6. The finite volume ef- 
fects can be seen from the resulting estimates of the scale: 
1/a = 1.87(6), 1.88(4) and 2.01(2) GeV, for 18^, 24^ and 
32'^ lattices, respectively. Thus we may have underesti- 
mated the scale by « 7%. It turns out, however, that this 
uncertainty is smaller than the range of scales resulting 
from the use of different physical quantities, and so can 
be subsumed into the quenching error discussed below. 

Data for the ai meson has much larger errors and we 
use a simple linear fit. The result, shown in Fig.|Bl gives 



aMa, = 0.58(4) -I- 2.1(1.0)(amg) 



(4) 



Again, the chirally extrapolated value is used to deter- 
mine the estimate of 1 /a given in Table |nj 

One of the well-known uncertainties introduced by 
quenching is that different physical quantities lead to dif- 
ferent values of the lattice spacing. Since we are inter- 
ested here in comparing with other quenched results for 
rUs and Bk, we follow most previous calculations and 
determine our central value for the scale, 1/a = 1.95 
GeV, using p masses. This lies within the range of values 
quoted above from the JLQCD Bk calculation, and is 
close to the value, 1/a = 1.855(38) GeV, they use when 
estimating rUs [13 • Nevertheless, to understand the im- 
pact of the scale uncertainty, and, as noted above, to in- 
clude possible finite volume errors, we also analyze subse- 
quent data using 1/a = 2.1 GeV. This value is consistent 
with our estimate from ai as well as the result obtained 
using the Sommer parameter rg IJH (1/a = 2.12 GeV). 
The latter is derived from the static qq potential and thus 
is independent of the fermion action. 



III. STRANGE QUARK MASS 

Our results for the masses of the lattice psuedo- 
Goldstone pion (spin-taste 75 ^5) are presented in Ta- 
ble lUll The results with different sources and sinks are 
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FIG. 1: EfTective mass plot of aAIp at quark mass 0.01. 
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FIG. 2: Effective mass plot of aAIp at quark mass 0.02. 



consistent, and we use the weighted average of the four 
results in the subsequent analysis. 

The strange quark mass rUs is determined by requiring 
a fictitious ss pseudoscalar mass to match the physical 
value of (2M|. — M^) which corresponds to {aMpsY = 
0.1234 with 1/a = 1.95 GeV. A linear fit for {aMpsf 
versus aruq works well, as shown in Fig. [3 This fit gives 
arris = 0.0520(7), and thus m, = 102(1.3) MeV. Re- 
peating the analysis with 1/a = 2.1 GeV leads instead 
to 110(1.5) MeV. We stress that these results are not 
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FIG. 3: Effective mass plot of aMp at quark mass 0.03. 
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FIG. 4: Effective mass plot of aMp at quark mass 0.04. 



TABLE III: Pion masses using axial and pseudoscalar oper- 
ators from Left {t — 10) and Right {t — 36) wall sources. 





aM.^{A4.,L) 


aM^{AA,R) 


aM^{P,L) 


aAU{P,R) 


mi 


0.1697(29) 


0.1682(30) 


0.1658(50) 


0.1644(56) 


1712 


0.2266(27) 


0.2255(28) 


0.2248(40) 


0.2224(42) 


ms 


0.2732(27) 


0.2725(26) 


0.2716(35) 


0.2695(34) 


7714 


0.3136(27) 


0.3134(24) 


0.3120(32) 


0.3106(30) 



very sensitive to the chiral fit form used. For example, 
a fit that includes a quenched chiral logarithm Il7| 
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FIG. 5: aMp vs. quark mass. 
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FIG. 6: aMai vs. quark mass. 



and is forced to pass through the origin reduces nis 
by 2.1(6) MeV. Such consistency is not surprising since 
arus — 0.052 is larger than the simulated points, whereas 
quenched chiral logarithms arc important only at masses 
smaller than aniq — 0.01. 

As discussed in the previous section, we expect that fi- 
nite volume effects should be small in our determination 
of TOs, because our heaviest two quark masses dominate 
the determination, and these have relatively large val- 
ues of Mt^L, 4.3 and 5.0 respectively. According to the 
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FIG. 7: (aM^)^ vs. quark mass. 



FIG. 8: Bk at quark mass 0.01. 



quenched chiral perturbation theory analysis of Ref . ^3 , 
one would expect to be larger than its infinite vol- 
ume value by 1-2% in this quark mass range. This would 
lead to our finite volume result for rus being 1-2% smaller 
than the infinite volume value. This estimate assumes the 
quenched hairpin parameter to be <5 « 0.2; using more re- 
cent values of (5 « 0.1 reduces the efi^ect proportionally. 
That the finite size effects are no larger than this size 
is supported by the numerical analysis of finite volume 
effects given in Ref. [T8| . 

Finite volume errors also enter into our result for nis 
through their effect on the scale, as discussed in the pre- 
vious section. We choose, however, to quote a value for 
nis for a definite choice of scale, so as to allow more 
straightforward comparison with other results. In par- 
ticular, using the one-loop matching factor from Ref. 
and the scale 1/a = 1.95 GeV, we find the renormal- 
ized mass m,(MS, 2 GeV) = 101.2 ± 1.3 ± 4 MeV. Here 
the first error is statistical, while the second is from the 
systematic effects that we control aside from the scale 
uncertainty. It is dominated by the uncertainty in Z,„, 
which we estimate as 4% by assuming a two-loop term of 
size ±1 X {as)^- It also contains the uncertainty from the 
form of chiral fit used, and from the finite volume errors 
in discussed in the previous paragraph. Note that we 
take the central value from the fit form without quenched 
chiral logarithms so as to better compare to the results 
of Ref. 0- 

This result for nis allows us to study the efficacy of 
HYP improved staggered fcrmions. The state-of-the-art 
quenched estimate for unimproved staggered quarks (ob- 
tained using the same definition of rus, and Afp(i) for 
setting the scale) is ms(MS,2GeV) = 106.0 ± 7.1 MeV, 
after extrapolation to the continuum limit . Our first 



observation is that our result at /? = 6.0 agrees with this 
continuum value, consistent with our expectation that 
errors should not be large. In this respect, we note that 
it was necessary to go down to lattice spacing a — 0.06 fm 
with unimproved staggered fermions in order to obtain 
the continuum estimate |lClj |. 

It is also useful to compare with the results from 
Ref. [l3j obtained at our coupling, /3 = 6. Their bare 
quark mass, arug = 0.0244 or rn^ ~ 45 MeV, is much 
smaller than ours. Using non-perturbative renormaliza- 
tion they find ms{MS, 2GeV) = 114 MeV. The very large 
matching factor, ~ 2.5, shows the need for non- 
perturbative renormalization with unimproved staggered 
fermions. Indeed, using one-loop matching they find the 
significantly smaller value TOs(MS, 2GeV) = 84 MeV. By 
contrast, our matching factor is very close to unity, illus- 
trating one of the advantages of HYP smeared staggered 
fermions. 

Quantifying improvement in discretization errors is 
more difficult. The unimproved (but non-perturbatively 
renormalized) result drops by 8 MeV between /? = 6 and 
the continuum, whereas our result is 5 MeV lower than 
the continuum value. Since these differences are compa- 
rable to the errors, the only definite conclusion we can 
draw is that the discretization errors appear to not be 
worsened by improvement.'^ 



This conclusion is bolstered by the fact that Ref. llCl uses a 
smaller scale than us (1.85 rather than 1.95 GeV). Had they used 



the larger scale, their final result at /3 
more from the continuum value. 



6 would have differed 
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FIG. 9: Bk at quark mass 0.02. 
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FIG. 11: Bk at quark mass 0.04. 
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FIG. 10: Bk at quark mass 0.03. 



IV. Bk 

The ratio of correlators corresponding to Eq.^is mea- 
sured on the interval 1 < t < 25 between two random 
U{1) sources ^19] placed a.t t — and 26. We find 
that the individual pseudoscalar meson correlators ex- 
hibit contamination from excited states up to f» 9 time 
slices from the sources. For this reason we choose to 
make constant fits to the central part of the plateau on 
time slices 9 < t < 16 even though the estimate is sta- 



TABLE IV: Results for bare Bk, BkCNDR,!! = l/a,L) and 
estimate of finite volume shift: 5Bk(L) — Bk{L = oo) — 
Bk{L). Errors are statistical. 



rriq 


bare Bk 


BK{l/a,L) 


3BK{l/a,L) 


nil 


0.514(61) 


0.416(57) 


+0.0094(38) 


7712 


0.614(26) 


0.542(25) 


+0.0013(5) 


m-j 


0.658(16) 


0.596(15) 


-0.0003(1) 


m4 


0.686(11) 


0.628(11) 


-0.0006(2) 



ble over the range 4 < i < 20. These fits are shown in 
Figs. ISIlTl In the first column of Tabic Hvl we give the 
resulting bare values for Bx, i-e- with all renormalization 
constants set to unity, for each of the quark masses. In 
the second column we give the results after renormaliza- 
tion to MS,NDR scheme at scale ji — l/a. 

To quote results in the MS scheme we use the one 
loop renormalization factors of Ref. , with the match- 
ing scale chosen to be g* = 1/a. The coupling as{q* = 
1/a) = 0.192 is calculated from the plaquette expectation 
value (P = 0.59367) using the method of Ref. At the 
physical kaon mass, which corresponds to arUq = 0.026, 
the one-loop corrections lead to a ~ 10% change in Bk- 
This is very similar to the corresponding shift with unim- 
proved staggered fermions. 

To extract Bk at the physical kaon mass we fit the data 
to the form predicted by quenched chiral perturbation 
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theory including finite volume corrections 



Bk = 6n< 1-6.0 



M 



K 



log 



Ml 
(4^/)2 



+ — [-25i(Mi,0,L) + M|.g2(Mi,0,L) 



(5) 



where / is the decay constant, which we fix to 132 MeV. 
The finite volume dependence enters through the func- 
tions gi, defined in Rcf. _21]. This dependence of 
B/f(NDR, /i — l/a) on the quark mass is shown by the 
solid line in Fig. E| with parameter values 6o — 0.23(9), 
bi = 0.5(1. l)/GeV2 and 62 = -1.2(1.3)/GeV''. Since the 
prediction for the finite volume corrections becomes un- 
reliable once MkL becomes small, we do not display the 
fit function below MkL = 2. 

Our results are consistent with the curvature predicted 
by the chiral logarithm in Eq. |SJl. Indeed, we can set 
62 = and obtain a good fit [with bo — 0.283(29) and 
bi = —0.30(19) GeV^^], showing that the curvature can 
be accounted for by the logarithm alone. Another consis- 
tency check is that bi and 62 agree with the expectations 
of naive dimensional analysis, namely ~ I62I ~ 1 
in units of the scale, ~ 1 GeV, of chiral perturbation 
theory. Taken as a whole, previous work is inconclu- 
sive concerning the presence of the chiral logarithm with 
predicted coefficient, largely due to the relatively high 
quark masses used (m > TOs/2). It is only by extending 



the range to 771^/5 that we find evidence, albeit not con- 
clusive, for the onset of the expected chiral logarithm at 
small quark masses. 

It is important to obtain a good fit to the chiral behav- 
ior in order to reliably extract estimates of finite volume 
corrections. The smallest value of MkL is 2.64, so one 
expects such corrections to be large ■ For Bk, how- 
ever, the cancellation between gi and 32 terms suppresses 
these corrections. Based on our chiral fit, the third col- 
umn in Table IIVI gives estimates for this finite volume 
shift. Note that the correction is non-monotonic in Mk, 
due to the cancellation noted above. The infinite volume 
prediction is shown in Fig. 1121 From this we conclude 
that for physical Mk the finite volume shift in Bk is 
much smaller than quoted errors even on our small lat- 
tices. This conclusion is supported by the absence of 
finite volume errors in the JLQCD results a,t P — 6 using 
unimproved staggered fermions [Tll | . 

Our final results are obtained by evolving from l/a to 
/i = 2 GeV using two-loop renormalization group running 
for iV/ [23. We find 



S/f(NDR,2GeV) 



Bk 

vRGI 
"0 



0.578 ±0.018 ± 0.042, (6) 
0.806 ±0.025 ±0.058, (7) 
0.314 ±0.124 ±0.176. (8) 



where BfP^ is the renormalization group invariant 
B— parameter [23, with b^^^ its value in the chiral limit. 
The first error combines that from statistics and those 
due to the chiral interpolation (or extrapolation for 69 )■ 
The second is our estimate of the uncertainty from us- 
ing one-loop matching factors explained below. Aside 
from the errors due to quenching and the use of de- 
generate quarks, which we do not address here, other 
systematics lead to changes smaller than the perturba- 
tive error. For example, using the scale from tq reduces 
BK(NDR,2GeV) and b^^^ by 0.025 and 0.007, respec- 
tively, while setting 62 = in the chiral fit increases them 
by 0.011 and 0.079. If we use f ^ Jk ^ 159.8 McV in- 
stead of / = 132MeV in Eq. El and fit the data, Bk 
changes by less than 0.01%, while 60 increases by 5%. 

The error associated with unknown corrections is 
estimated as follows. We can write Bk — Ba+Bv, where 
V and A refer to vector-vector and axial-axial parts of 
the operator in Eq. Ba,v can each be decomposed 
into one and two color-trace parts |j] . Each of these four 
components of Bk is proportional to log(M|-) and thus 
diverges in the chiral limit, although their sum does not. 
Using one-loop matching there is an incomplete cancel- 
lation, and the resulting Bk should diverge in the chiral 
limit, although this feature is expected to manifest itself 
at much smaller quark masses than studied here. Indeed, 
the bare values for By, a do indicate a divergent behav- 
ior. Because of the residual divergence, we cannot simply 
estimate the error, in particular in 60, by multiplying by 
an overall relative correction of ±a{q*)'^ (as we did for 
iris). Instead, we recalculate Bk after adding zLa{q*)^ to 
the matching factors for each of the four components of 
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TABLE V: Bk(NDR,2 GeV) at /3 = 6 for gauge invariant 
(GI) and non-invariant (NGI) operators, before and after re- 
moving the fitted c? and c? terms. Data from Ref. |ll|. 



Type 


Uncorrected 


removed 


c? removed 


GI 


0.6790(16) 


0.55(7) 


0.76(7) 


NGI 


0.7128(14) 


0.61(7) 


0.73(7) 



Bk in turn, and take the largest variation as the error. 
The resulting uncertainty, quoted above, is larger than 
the statistical error and, as expected, grows rapidly in 
the chiral limit. 

Even though we need more high precision data at 
lighter quark masses to pin down the chiral extrapo- 
lation, it is nevertheless interesting that our estimate 
h^Gi ^ 0.314 ± 0.124 ± 0.176 is in good agreement with 
recent estimates 0.29(15) and 0.36(15) ^ obtained 
using 1/A'c expansion. 

We now compare our estimate with the state-of-the- 
art results obtained by the JLQCD collaboration 
using unimproved staggered fermions and argue that 
HYP smearing reduces discretization errors. The JLQCD 
result in the continuum limit is (NDR, 2 GeV) = 
0.628 ± 0.042. Our first observation is that our result 
at = 6 is consistent with this continuum result. On 
the one hand, this indicates that the c? errors with HYP 
fermions are not large, as for nis. On the other hand, 
the value of Bk at /? = 6 with unimproved staggered 
fermions [0.679(2)] is also consistent with the continuum 
result, giving no evidence of improvement. 

We can go further, however, using the details of the 
continuum extrapolation provided by Ref. |lll |. Their fit 
included both and terms. Using their fit param- 
eters we can determine two additional estimates of Bk 
at /3 = 6.0: removing only the Oio?) term (and not the 
O(a^) discretization correction) and vice- versa. The orig- 
inal results and those corrected for either the discretiza- 
tion or perturbative errors alone are given in Table [V] 
The point we wish to make is that in the case of gauge 
invariant operators (which are those we use) the JLQCD 
fits imply that, at /? = 6, the total result 0.68 contains an 
O(a^) contribution of ^ 0.13 and an 0(o?) contribution 
of ~ —0.08. (Corrections for non gauge-invariant oper- 
ators, which are not the operators of choice, are some- 
what smaller but show the same pattern.) Thus, in a 
formulation where only discretization errors were elimi- 
nated or substantially reduced one should expect a final 
result closer to 0.55 at a w 0.1 fermi. Our estimate with 
HYP smearing, Bk(NDR,2 GeV) = 0.58(4), is indeed 
consistent with this and significantly different from the 
unimproved JLQCD result 0.679(2). This suggests that 
discretization errors have been reduced by using HYP 
smeared fermions. 



This conclusion is supported by the fact that all cal- 
culations with domain wall or overlap fermions, which 
are also expected to have small discretization errors and 
small perturbative corrections, find values at /? = 6 
consistent with ours: 0.575(6) (Ref. jH), 0.532(11) 
(Ref.JgTf). 0.563(31) (mean of data at /3 = 5.9 and 6.1 in 
Ref. lUl), and 0.63(6) (Ref. '29ll|). Here, only statistical 
errors have been quoted. 

Finally, we consider the results in TablelVlwith "a^ re- 
moved" . These correspond approximately to using non- 
perturbative matching factors, and should thus expose 
the "true" a} errors in the unimproved results. Unfor- 
tunately, the large errors preclude definitive conclusions. 
Nevertheless, the fact that our result 0.58(4) differs from 
the "a^ removed" unimproved result of 0.76(6) by about 
2(7, while lying closer to the continuum result 0.62(4), is 
consistent with our conclusion that discretization errors 
are reduced by using smeared links. 



V. CONCLUSION 

We conclude that improved staggered fermions are a 
viable and promising option for calculations of TOs and 
Bk in full QCD simulations. The difficulties observed 
with unimproved staggered fermions (ill behaved pertur- 
bation theory for Z„j in the case of TOs and discretization 
errors in Bk) are greatly reduced. Our study suggests 
that reliable calculations should be possible on the en- 
sembles of lattices being generated with dynamical im- 
proved staggered fermions without requiring very small 
lattice spacings. There are two caveats, however. To 
reduce the uncertainty due to the two-loop term in the 
renormalization constants below our estimates of 4% in 
TOs and 7% in Bk will require a demanding two-loop or 
non-perturbative calculation of matching factors and the 
calculation of a larger set of lattice matrix elements. Sec- 
ond, in the case of Bk, the estimate in the chiral limit is 
very sensitive to errors in the matching factors, as well 
as to the chiral extrapolation. 



VI. ACKNOWLEDGEMENT 

This calculation has been done on Columbia QCDSP 
supercomputer. We thank N. Christ, C. Jung, C. Kim, 
G. Liu, R. Mawhinney and L. Wu for their support on 
this staggered e'/e project. This work is supported in 
part by BK21, by Interdisciplinary Research Grant of 
Seoul National University and by KOSEF contract ROl- 
2003-000-10229-0, the US-DOE grant KA-04-01010-E161 
and contract DE-FG02-96ER40956. We gratefully ac- 
knowledge discussions with A. Soni. 



[1] D. Becirevic, Nucl. Phys. B (Proc. Suppl.) 129 & 130, [2] C. T. H. Davies et ai, Phys. Rev. Lett. 92, 022001 (2004). 
34 (2004). 



9 



[3] G. Kilcup, R. Gupta and S. R. Sharpe, Phys. Rev. D 57, 
1654 (1998). 

[4] M. Golterman, Nucl. Phys. B Proc. Suppl. 73, 906 
(1999). 

[5] J. F. Lagae and D. K. Sinclair, Phys. Rev. D 59, 014511 

(1999); G. P. Lepage, ibid. 074502. 
[6] W. Lee and S. Sharpe, Phys. Rev. D 66, 114501 (2002). 
[7] A. Hasenfratz and F. Knechtli, Phys. Rev. D 64, 034504 

(2001). 

[8] K. Orginos and D. Toussaint [MILC collaboration], 

Phys. Rev. D 59, 014501 (1999). 
[9] W. Lee and S. Sharpe, Phys. Rev. D 68, 054510 (2003). 
[10] S. Aoki, et ai, Phys. Rev. Lett. 82, 4392 (1999). 
[11] S. Aoki, et al, Phys. Rev. Lett. 80, 5271 (1998). 
[12] C. Sui, PhD. Thesis, Columbia University (2000). 
[13] M.F. Golterman, Nucl. Phys. B 273, 663 (1986). 
[14] M. Booth, G. Chiladze and A. F. Falk, Phys. Rev. D 55, 

3092 (1997). 
[15] R. Sommer, Nucl. Phys. B411, 839 (1994). 
[16] C. W. Bernard and M. F. L. Goherman, Phys. Rev. D 

46, 853 (1992). 



[17] S. Sharpe, Phys. Rev. D 46, 3146 (1992). 
[18] S. Aoki, et al, Phys. Rev. D 50, 486 (1994). 
[19] D. Pekurovsky and G. Kilcup, Phys. Rev. D 64, 074502 
(2001). 

[20] C.T.H. Davies, et al, Phys. Rev. D 56, 2755 (1997). 
[21] J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987); 

188, 477 (1987). 
[22] D. Becirevic and G. Villadoro, Phys. Rev. D 69, 054010 

(2004). 

[23] J. Buras, M. Jamin, and P. Weisz, Nucl. Phys. B 347, 
491 (1990). 

[24] J. Prades, J. Bijnens and E. Gamiz 

arXiv:hep-ph/0501177 
[25] O. Cata and S. Peris, JHEP 0303, 060 (2003). 
[26] A. Ah Khan, et al, Phys. Rev. D 64, 114506 (2001). 
[27] T. Blum, et al [RBC Collaboration], Phys. Rev. D 68, 

114506 (2003). 
[28] Thomas DeGrand, Phys. Rev. D 69, 014504 (2004). 
[29] Nicolas Garron, et al, Phys. Rev. Lett. 92, 042001 

(2004). 



